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ABSTRACT 



Conjugate basis functions which are continuous through- 
out a model are used to approximate a stress field. This 
conjugate approximation has less mean error than approxima- 
tions calculated using conventional finite element methods 
and is a better approximation at extreme values of stress. 
This thesis consists of two major parts. The appropriate 
analytic expressions for conjugate stress fields are set 
forth and applied to the constant strain triangle finite 
element. Analyses are then performed with these new stress 
calculations and the results compared to some previous 
stress calculations. 
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INTRODUCTION 



I. 



A. BACKGROUND 

In a recent paper Oden and Brauchli [Ref, 1] state that 
one of the difficulties encountered in using the finite 
element method to approximate the stress field in a model 
has been that the stresses calculated are at best only 
averages over each element and, in general, are discontinu- 
ous between elements. Continuous stresses can be obtained 
only if strain degrees of freedom (DOF) are introduced at 
each node. This is achieved at great computational expense. 
Some of the schemes which are used to calculate the state 
of stress at a node are mentioned in Ref. 1 and another 
method which is based on conjugate approximations is pro- 
posed. In Ref, 2 Oden expands on this method explaining 
the theory of modeling a function in general and using con- 
jugate space functions for modeling in particular. 

Briefly the methods of evaluating nodal stresses re- 
ferred to in Ref. 1 include taking an average of the 
stresses of each element incident at the node; using a 
weighted average method developed by Wilson [Ref. 3]; using 
averages based on equating nodal forces and element stresses 
in a procedure developed by Turner, Martin and Weikel 
[Ref, 4]; or using stiffness matrices proposed by Gallagher 
[Ref. 5] which are developed relating both the stress and 
displacement fields of each element. Among all methods, the 
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proposed method, using conjugate bases, has been shown to 
give the best fit in a least squares sense. 

B. SCOPE AND OBJECTIVES 

In this thesis the conjugate approximation theory is set 
forth in detail and is applied to illustrative examples. 

The appropriate analytic expressions for conjugate stress 
fields are derived for the constant strain tiangle (CST) 
finite element and incorporated into a CST computer program. 

Numerous examples of in-plane problems are worked com- 
paring both the CST solution and the conjugate basis method 
(CBM) solutions and the exact solution if it is available. 

Finally, some conclusions are reached concerning the 
type of in-plane problems where the CBM may be used success- 
fully and what the cost in additional computer time may be 
for this continuous best solution. 
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II. CONJUGATE APPROXIMATION THEORY 



A. MODELING A FUNCTION 



A function F(x) over a region or domain R can be 
modeled by another function F(x) in the domain R where 
F ~ F and R z R. The region R is subdivided into elements 



where each element has two or more nodes depending on the 

degree and type of the approximation function over the 

element. Within R all elements have the same number of 

nodes and each node in the global model R is designated 

X^, (A = 1,G), where G is the total number of system nodes. 

The domain of the element e is r and each node within this 

e 

N 

local model is designated x^, (N=l, N^) , where N^ is the 

number of nodes in an element. 

For example, a one-dimensional axial bar may have the 
following three element representation (here R = R) . 





2 



x 



e 



L 



Figure 1. Modeling a Bar. 
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Over each element the function F (x) is represented by a 

0 

local function, f (x) , which is only defined on the element 

N N N N 

e. Its value at x is denoted f , that is, f = f (x ). 

e e e e 

Hereafter, the index e denotes element e. In the figure 

N 

below, the functions F(x), and f^ for element 2 (i.e. e=2) 

are shown. 





Figure 2. Global and Local Functions. 

B. INTERPOLATION FUNCTIONS 

A set of independent local interpolation functions 
i^(x), (N=l, N^) , is defined over each element such that 

0 

i|^(x) = 1 at node N and zero at all other nodes. (For ele- 
ments with linear displacement fields, such as the elements 
used in this thesis, is linear between nodes.) For a 

one-dimensional 2 node element, a set of interpolation func 
tions for representation of a linear field is, 

t \>1 = i - e ; <J>2 = c (1) 

where £ is a local coordinate over the closed domain [0,1] 
and origin at local node point 1. 
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Figure 3. Local Interpolation Functions 

0 

The local function, f (x) , is approximated over the ele- 

— 0 

ment by f (x) where 

f e (x) = N = l, N g . (2) 

— e N e 

[The convention f (x) = f ^ (x) means summation over 
repeated indices, that is 

f 6 (x) = + f ^2 ( x ) + • • • + ^ Cx) • 

e 

This summation convention will be used throughout,] 

_ e N 

[Note that at a node the value of f (x) is f^.] Here 

2 -2 

the local function f (x) is modeled by f (x) with linear 

2 2 

interpolation functions, ^ and 
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Figure 4. Approximation to Local Function 



C. INCIDENCE OPERATORS AND ASSEMBLY 

e A e N 

The function and ft^ are known as incidence operators 

and are used in relating the global model R to local ele- 



ments r . 

e 


They have the value of 


1 or 0 as 


follows 




i 


if node A in R is coincident 


with 


e A 




node N in r in the 


ass emb led 


model. 


» A 

a n = 




e 






0 


otherwise 







(3a) 



e 

ft 



N 

A 



1 if node N of element r is coincident 

£ 

with node A of the assembled model, (3b) 

0 otherwise 



A is called a mapping function because it gives the value 



of a local node to the incident global node. Likewise, ft 



is called a decomposition function because it gives the 



global value to 

f a = > f N 

N e 



f N - !)■? F A 
e A 



the corresponding local node. 

(Note this is summed over both 
e and N . ) 



(4) 

(5) 
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The model of a function F(x) in terms of 4^ element in- 



terpolation function can now be given by 

F (x) ~ F(x) = ^F A E=l, N 



( 6 ) 



e = l 



N = l, N 

< 

A=l, G 



where N is the number of elements, 
r 

Equation (6) is true everywhere except at the global 

nodes with m elements coincident. Here the value of F(x) 

is mF(x) times greater than it should be. Accordingly, at 

the nodes Equation (6) becomes 

E_ 

e 



5(x) - i E *n' x) s a fA 



(7) 



e = l 

We note that Equation (7) gives continuity of F at 
nodal points . 

A set of global interpolation functions 4>^(x) are now 



defined as follows 

<}> A (x) 



E *N <X) “a 

Te=l 



* e ♦: 



(x) ft 



N 



except at nodes 



at nodes (8) 



e = l 



It will be demonstrated in the examples that <J)^(x) is 

a linear combination of the element interpolation functions 

e e 

4^ 1 s . Moreover, since the local 4^ are independent, the 

global 4^ functions are independent also. There is a unique 

4^ for each global node. If the element interpolation 

functions are independent, the 4^* s are independent func- 



tions. Now 



F(x) = c|> A (x) F* - ^ ^n ( x) f e 

e = l 



(9) 
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The domain of a global interpolation function <j)^ extends 
over all elements incident at X^ . This function has the 
value 1 at and zero at all other nodal points. For ex- 
ample, with the previous linear interpolation functions, 
and the three element axial domain, we have 




Figure 5. A Global Interpolation Function. 

D . FUNDAMENTAL MATRIX 

The set of G linearly independent functions { } form 
the basis for a G-dimensional space $ where G is the total 
number of global nodes. The fundamental matrix of the space 
$ is defined as 

c A r E A> r = 1 , G (10) 

The notation use< ^ to denote the inner pro- 

duct of (J>^ and p . Here the inner product is taken as 

«J> A ,<t> r > - / R <f> A 4> r dR A,r = 1 , G (11) 
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is a positive, definite 



Since / R (f> A (J> A dR > 0, C Af 
matrix and its inverse can always be formed. It is defined 
as 

-1 Ar 

(C Ar ) = C (12) 



This will be the transformation matrix used to form a 
new set of basis functions in the G-dimens ional space. More- 
over, since f <f) ^ <(> p dR = cf) p (J)^ dR, C^p is also a symmet- 
ric matrix. 

It is convenient to construct C^p from local interpola- 
tion functions as follows. Substitute Equation (8) in 
Equation (10) to obtain 

E E 



e, T f , 



c Ar ^a’V ^ ^ Q r < ^n» V 



(13) 



e=l f=l 



e f 

Here is the inner product of the local interpola- 

tion functions. It is non zero only if e=f since each func- 
tion is only defined on its own element. The inner product 
e e 

<} Pxi> > is called the local fundamental matrix for element 
N m 

e 

e > C NM * 



Equation (13) now becomes 

' E 



'Af 



•E 



e M 6 
c nm 



(14) 



e=l 



E. CONJUGATE SPACES 

A conjugate basis is one which is reciprocally related 
and interchangeable as to properties with the original 
basis. A set of G linearly independent functions on R, 
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{<J) A .}, a so-called conjugate basis set, must satisfy the re- 
ciprocity relation, 

where 6. A = A ^p (15) 

This means that the sets {(J)^} and form a countable 

biorthogonal basis of the space Any sets and {<j)^} 

which satisfy Equation (15) are referred to as conjugate 
basis of $. 

In order to define {(f)^} uniquely, they are formed as 
follows 

4> a (x) = c Ar 4> r (x) (16) 

It can be shown that the (f)^*s are independent functions, 



and that 



Moreover 




> 4> r > 



j kt 



qy A (x) = C Ar <J> r (x) (17) 

< c AK <f> K * (f r > = c AK <<$> k a t > 

<Sp as required by Equation (15) • 



Using the incident operators defined in Equation (4), 

the local conjugate basis interpolation functions are 

,N .A /10 v 

± e = ^ (18) 

and making the substitutions from Equations (16) and (9) 

gives 



. N _ n AT . 

~ ^ A ^ 

n Ar v , f, n £m 

fi, C S ip (x) fip 

f = l 



(19) 



A,r = i , g 



H = 1 , N (20) 



Equation (20) highlights the most important feature of the 
local conjugate basis functions. The functions for each 
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node of each element are linear combinations of the basis 
functions of all G nodes. This gives each local node a 
basis which is continuous over the entire model. The first 
example gives a numerical illustration of this. 



F, APPLICATION OF THEORY TO FINITE ELEMENT 

E 

Recall Equation (10) where F(x) = F^<J>^= E ^ f^ • An 

e = l p 

alternate representation of F(x) in terms of the <J> basis 
is obtained as follows 

F = f\ = f a c Ak 4> K 



_ , K v , N _e 

= V ■ te f N 
e = l 



( 21 ) 



Here F is the value of F(x) with respect to the con- 
K 

jugate basis. It is calculated by forming the inner product 

of F and <J> . From Equation (21) we have 

K 

<F -V ■ <F /‘-V ■ 



= F 

N K 



= F 



K 



( 22 ) 



The local values of F are obtained by the usual decom- 
position, 

e A e A _ _ „ 

(23) 



£ n - «N' F a - «N <F -V ■ <f -*n > 

f N = ft A f a = fi A <F,4) A > = < F , jjj N > 
e N N * Y ,x e 



(24) 



Now the original function can be approximated in the 
conjugate basis using Equation (21). This is the procedure 
used in solving problems using the conjugate basis method 
(CBM) . 
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G. BEST APPROXIMATION 

Before proceeding, it is important to verify that ap- 
proximating a function using the conjugate basis functions 
will give the best approximation in a least square sense. 

This is proven simply and clearly with a simple example. 

2 

The function F(x) = x is to be approximated in the in- 
terval x = 0 to x=l by a linear polynomial F(x). Here bases 
of F are functions $^ = 1 and (f^-x 



F = a^<f>^ + a 2 <}> 2 = a ^ ( 1) + a 2 (x) (25) 

The difference between the actual function and its 
approximation is called the error, E 



E = | F - F | 



(26) 



The function J = J(a.) is a measure of the error over 

1 

the interval. 

J(a.) = E 2 dx = /"'’[(a, + a 9 x) - x 2 ] 2 dx 
1 o o 1 z 

t / \ 2 , 2 ,121 ,1 

J(a.) = a;L + ai a 2 ~ ^ a l + 3 a 2 ~ 2 a 2 + 5 (27) 

The error is minimized by setting the partial deriva- 
tion of J (a^) , with respect to a_^, equal to zero and solv- 
ing the resulting set of equations for a_^ f s. 



Now a^ 



1 

6 ; 



3J(a 1 ) 
3 a, 



= 0 = 



2a l + a 2 - 3 



8J(a 2 ) 

9a» 



= 0 = =1 + |a 2 - i 



(28) 



a „ = 1 



and F = - t + x (29) 
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Using the conjugate basis method with basis functions Jp^ 9 
F is determined as follows: 

1. Assume the local basis functions are ^ = 1 - x 

and 1^2 = x. 

2. Form the local fundamental matrix c XT = 

NM N M 

Here we obtain, 



C NM 



— [ 2 1 ] 
6 l l 2 J 



(30) 



3. Take the inverse 



NM .4-2, 
c ' - 2 4 



(31) 



N NAT 

4. Form the conjugate basis function Jp^ = C <p ^ 

in this example <})p = because there is just one element 

.1 ^11 . , 0 12 . . 
ip = C (p^ + C ^2 = 4 - 6x 



2 21 22 
Jp = C <p ± + C <P 2 = - 2 + 6x 



(32) 



e e 

5. Determine a. where a. = <F, \p . > 

l i r i 



a 1 = f 1 x 2 (^ 1 dx = j-j; a 2 = f 1 x 2 (p^dx = ~ 

o o 

1 e ^ 

6. Approximate F(x) with F = E a. \p i = l,2 

e=l 1 6 

F = ^-(4-6x) + |(-2+6x) 



= " 6 + x 



(33) 
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Figure 6* Best Approximation 

This is the same approximate function obtained using the 
previous calculus method. Therefore, the best approximation 
of a given function F (x) in the subspace <j) is the fujiction 
F(x) whose components are the inner product of F(x) and the 
conjugate basis functions. Reference 2 contains a formal 
proof that the CBM is the best fit in the least square sense. 
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III. APPLICATIONS TO DEMONSTRATIVE EXAMPLES 



In all of the problems considered in this thesis, the 
displacement must be in the domain of the structure. For 
one-dimensional problems, the model is a bar displaced along 
the x axis. The two-dimensional model is a plane with dis- 
placements in either the x or y direction or both. 



3 



x 



y 



X 



Figure 7. One- and Two-Dimensional Domains. 

Appendix A contains a brief summary of the stress- 
strain and strain-displacement relations for the in-plane 
problem . 

A. ONE-DIMENSIONAL PROBLEMS 

1 . Nonlinear Displacement and Linear Properties 

To demonstrate certain properties of conjugate ap- 
proximation functions, example 9.2 in Ref. 2 will be 
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explored in detail. 

a. Statement of Problem 

Determine the stress in a nonhomogeneous , 3 unit 
length bar based on approximate displacement fields. The 
stress is given by 



o(x) = MX) 

where k(x) = k^Cl+x) 

2 

is the material modulus in units of force/length 

is the displacement field. Assume that 

2 

u(x) = a [ 1 - ] 



(34) 

(35) 
and u(x) 



(36) 



where a is a small constant and x is measured in units of 
length . 



With this displacement field, 
distribution is 



2ak 

a(x) = - ■ 9 — x(l+x) 



the exact stress 
(37) 



The model R consists of three one-dimensional 

elements, each of unit length, that is R = R = [0,3] = 

3 

U ft 
i e 
e=l 



(i) 



( 2 ) 



(3) 



Figure 8. Model for Problem III A1 . 
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The local basis functions, for a linear field, 
corresponding to a typical element, e, are the interpolation 
functions shown with Equation (1) , 

= 1 - £ and = S 

b. Global Interpolation Functions 

The global interpolation functions are calcu- 



lated using Equation (9), 

3 e 



N 



<f» A (x) = E ft” ^N (X) 
e = l 



A = 1,4 
N = 1,2 



NOTE: The notation (x) will be omitted unless it is needed 

for clarity. 

The above expression gives, 



;i , 1 



\2 , 1 



;i , 2 



= 1 + o 



n 1 , 1 

0 

ii , i 



> 2*2 

■ *2 
;2 , i 



+ o 

2 



‘^l 

I 

;i , 2 



i , i 



+ o 

l 



2 , 1 



+ 0 
2 



1 , 2 



= 0 +0 +0 

These are shown in Figure 9. 



2 2 2 




3 1 3 




3 2 


^1*2 


+ 




+ 


Q ± ip 


0 


+ 


0 


+ 


0 


2 . „ 




3 




3~ 


2 2 
^2 * 2 


+ 




+ 


2 

^2* 


0 


+ 


0 


+ 


0 


2 2 2 




n 1 3 




3 2 


n 3 \p 2 


+ 




+ 


S* 


1 *2 


+ 


1 *1 


+ 


0 


2 o o 




3 1 3 




3 ? 


«4*2 


+ 


Vl 


+ 




0 


+ 


0 


+ 


i ip 



(38) 
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*1 
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Recall that using Equation (8) the value of <J>^ 
at a node where 2 or more elements are incident is still 1. 
c. Local Fundamental Matrix 

The local fundamental matrix (LFM) is formed 
using the local basis functions as discussed in Section II-D. 



e ^ . e i e 

c — = < ^N * V 



'NM 



N , M = 1,2 
e = 1,3 



cn = / 1 ^(0^(Od£ = / 1 (1-C) (l-£)d£= | 
o o 

*12 - l 1 U-C)£ « - I 



L, - = /J(?xi-£)d£= | 



21 



■1 , 1 



•22 



= = r(0(0 dc = i 09) 



M 1 if 2 i, 

NoW C NM = 6 [ 1 2 



(40) 



Since the elements are geometrically identical, we have 

1 2 3 

C NM C NM C NM* 

d. Global Fundamental Matrix 

The global fundamental matrix (GFM) is formed 

using Equation (14), 

C 



^ e N e M e 
tA „ - Z ft, ft" c 
AT . A T NM 
e = l 



N , M = 1,2 
A,r = 1,4 
E = 3 



Since C^p is a 4 x 4 symmetric matrix, a total 
of ten separate calculations would be needed to determine 
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all of the elements of the GFM. 



The calculations for c 



11 



and 



C 22 are s ^ own t> e l° w 



C 11 


= 


^1^1 


3 c-n+ 


^1^2 


"l2 + 


1 2 1 1 


^21 + 


1 2 1 2 

Si SI 


1 

C 2 2 




+ 


2 2 
11 


‘ll + 


2 2 
12 
Sl^Sl^ 


‘l2 + 


2 2 
2 1 


*21 + 


2 2 
2 2 
Sl^Sl^ 


2 

C 2 2 




+ 


3 3 
11 


^11+ 


3 3 
12 


i i2 + 


3 3 
2 1 


^21 + 


3 3 
2 2 
Sl*Sl* 


3 

c 2 2 




= 


£[i.i 


♦2 +0+0+0 


+0+0+0+0+0 


+0+0+0 ] 






= 


2 

6 
















C 2 2 


= 


^l^l 

«2 n 2 


=n + 


^1^2 


c 12+ 


1 2 1 l 

^2^2 


^21 + 


•^2^2 

^2 fi 2 


1 

C 22 




+ 


^ 1^ 1 

Sl\Sl^ 


I + 
11 


^1^2 

«2 fl 2 


^12+ 


^2^1 

*2*2 


^21 + 


2 2 
2 2 

^2^2 


2 

C 2 2 




+ 


^ 2 ^ 1 
°2 fl 2 


C ll + 


3 3 
J 2 
Sl\Sl^ 


C 12 + 


^2 3 1 


c 2i+ 


^2 3 2 

^2^2 


c 2 2 



below . 



= y -[ 0+0 + 0 +2 +2 + 0+0 + 0 + 0+0 + 0+0 ] 

o 

4 

6 

The complete GFM and its inverse are shown 



(41) 



'Ar 



i 

6 



2 1 0 0 




2 6 -7 2 -1~ 


14 10 


„Ar 6 


-7 14 -4 2 


0 14 1 


C ‘ 45 


2 -4 14 -7 


0 0 12 




-1 2 -7 26 



(42) 



In practice, the entire set of calculations 



using incidence operators is avoided and the LFM elements 



are placed in the GFM according to the following scheme. 

Each column and each row of the GFM is asso- 
ciated with one and only one global node. Likewise, each 
column and each row of the LFM is associated with just one 
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global node. In assembling the GFM, the ij element from the 
LFM associated with local nodes i,j is added to the k,l GFM 
element, where global nodes k,l correspond to local nodes 
i,j. If several elements are incident at a global node, 
then the appropriate elements from each of the connecting 
elements 1 LFM will be summed in the GFM, This procedure 
avoids numerous calculations where the incident operator is 
zero , 

e. Local Conjugate Basis Functions 

The local conjugate basis functions are calcu- 
lated using Equation (19), 

N e N A F c - 1 , 3 

£;(x) = C n <Mx) N = 1,2 

e n A,r = i,4 

The calculation for node 1 of element 1 is 

shown below, 

.1 Jl 11 , . Jl 12 , Jl 13, , *1 14, 

^-1 = c + c 4 > 2 + c <f > 3 + c $4 

h 21 , h 22 , h 23, ^ il 24, 

+ ^2 c ^ c ^2 + z C ^3 + ^2 C ^4 

ll 31,. . il 32, il 33, . il 34, 

* ub ^ e '1 ^3 c 9 2 3 c *3 3 c 9 4 

^1 41 42 ^1 43 ^1 44 

+ Q 4 c 4> ± + c 4 > 2 + c J (}) 3 + c (J ) 4 (43) 

For node 1 of element 1 ^ ^° r ^1 

0 for all other 

AT 1 

Therefore, substituting for C and (j>^ above becomes 

il = fj[26 4 ^ - 7(4>2 + 4^) + 2(4^ + 4^J) - 4 ^ 2 ] <**> 

where the \p ' s in brackets are the local basis functions de- 
fined in Equation (1). 
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Notice that only one row of the above calcula- 

e N 

tion was nonzero - that is the row where = 1. Using this 
fact reduces the actual calculations in a manner similar to 
the method used in forming the GFM. 

Consider the nonzero row of the above calcula- 
tion, Factoring out the common incidence operator leaves 
four terms which are the products of the V th element of the 
A row of the GFM inverse times (j)p • Recall that at a node 
<J)p = 1 and the local conjugate basis functions X = be- 

come equal to the A row of the GFM inverse. This greatly 
simplifies computer programming of the conjugate basis 
method . 

Returning to the example at hand. Global node 2 

is node 2 of element 1 and node 1 of element 2. Therefore, 

1 2 ^1 ^2 ^1 

^2 = ^2 = ^ an< ^ ^A = ^A ~ ^ f° r A = 1, 3, or 4. 

Showing only the nonzero terms of the calcula- 

t ion, 

tl - tl - c 21 4» 1 + c 22 4» 2 ♦ c 23 * 3 + c 2 \ 

- + 1«<4>2 + + + 2 *2 ] (45) 

A similar argument for global nodes 3 and 4 

gives 

±1 - ±1 = - 4(^2 + + 14(^2 + ( 46 ) 

and 

±1 = + 2(^2 + "7(^2 + + 26^2 3 (47) 

The local conjugate basis functions are shown 
graphically in Figure 10. Notice that each function is 
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52 







Figure 10, Local Conjugate Basis Functions. 
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continuous over the domain, x, and that the effects of nodes 

remote from the local node being considered diminish as the 

distance from that node increases. 

f. Conjugate Nodal Components and Global 
Representation 

Calculate the conjugate nodal components of the 
element stress field according to Equation (23). 

- <a(x) , 

0 

Since are in terms of 4> Cf(x) will be transformed to 

functions of 4 also, 

0(?) = + 24 2 for 0 < x < 1 

= 24 ^ + 64+4 for 1 x _< 2 

= 24^+104+12 for 2 x 3 

ak 

a factor of - — — has been factored out to enable compari- 
son with the results published in Ref. 2. Now 

o\ = f 1 U5 + 25 2 )(l-C)dC = | 

o 

o\ = f 1 (25 + 2^ 2 ) (5)d5 = | 

- o 

o 2 = Z 1 (2C 2 + 65 +4) (l-5)d5 = ±2- 

o 

a 2 = Z 1 ( 25 2 + 65+ 4) (5)d5 - f 2 

O 

o? = Z 1 (25 2 +i05+i2) (i-5)d5 = 

1 0 6 

o\ = Z 1 (25 2 + 105+12) (5)d5 = I 5 - (48) 

Now the global stress field can be modeled using Equation 

(21) > 3 e n 

a(x) = E a N £ e (x) N = 1,2 

e = l 
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At the global nodes we have 



0(0.) 



3 52 7 14 19 14 27 4 47 4 5_9 2_ 

6 15 6 15 6 15 6 15 6 15 6 15 



0 ( 1 . ) 



- 0.33 

3 (-14) 28, l 19.. 

6 15 15^6 6 ' 



8_ ( 27 47 4,59. 

15 ^ 6 6 ; 15 ^ 6 ; 



= 3.43 



( 49 ) 



Similar calculations are made for x=2 and x=3. 



The results are 



0(2.) = 11.7 
0(3.) = 23.6 

These results are plotted on Figure 12 and are 
also shown in Table 1, Since the exact stress field was 
known, a comparison of accuracy is made, 

g. Conventional Finite Element Solution (C.F.E.) 

Here we suppose the continuous displacement 
field U(x) is approximated by the piecewise linear field, 
U(x) where 

U(x) = ^-[ 9 4 > 1 (x) + 8^2 (x) + 5<f> 3 (x)]. (50) 



The 4> ' s are shown in Figure 9. This linear 
field was constructed to give U(x) = U(x) at the nodal 
points by evaluating U(x) at the nodes. The strain is now 



dU(x )_ 

dx 




d^ 

dx 



+ 8 




5 



d^ 

dx 



d<f > 2 

dx 

dcj>3 

dx 



1 



3 



5 



0 < x < 1 

1 £ x £ 2 

2 < x < 3 (51) 
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The conventional solution is obtained from 



a (x) 



k(x) 



dU (x) 
dx 



where the value of 



dU (x) 
dx 



is 



given above, 

2 . Linear Displacement Model and Linear Properties 

In the first example the exact stress function was 

0 

known and was used to calculate CL T (x). Consider the re- 

N 

suits if the stress field is approximated by 



a(x) = k (x) 



dU (x) 
dx 



(52) 



where ^ ^ X - is a constant over each element, 
dx 

This is the constant strain approximation. Since 

a constant strain triangle (CST) finite element computer 

e 

program is used to provide the T s in the two-dimensional 
problems considered later, using = constant provides a 
guide to the error which may be expected. 

U(x) as given in Equation (36) is calculated at 



each node and plotted vs x below assuming linear displace- 
ment between nodes. The slope of each line segment is 
and is shown in parentheses. 
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U(x) 




Figure 11* Assumed Linear Displacement Model 



The conjugate nodal components, 



(x) , are now 



given by 

Ot k 

a^(x) = - — ~ [ ( 1 +x) ij>®(x) dx (53) 

The results of the integration, disregarding the 
ak 

- — 5 — factor, are 



a 



l 

l 






a 



1 

2 




24 

6 




(54) 



The conjugate function representation of stress is 
given by Equation (21) and the results are plotted on 
Figure 12 and included in Table 1. 
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Figure 1,2. Exact Stress Distribution and Approximations 
to the Exact Stress. 
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TABLE I 



Summary of Results for Problems III A1 and 
III A2. 



a(x) 
x Exact 


a (x) 
CBM1 


Er ror 

% 


a (x) 
CBH2 


Error 

% 


C.F.E. 


a (x) 

C . Ave 


Error 

% 


0. 


0 


-0.33 


- 


-0.13 


- 


1 


i 


- 


1. 


4 


3.43 


14.25 


4.26 


6.5 


2,6 


4 


0 


2. 


12 


11.7 


2.50 


12.06 


0.5 


9,15 


12 


0 


3. 


24 


23.6 


1.66 


21.46 


10 .58 


20 


20 


16.67 



NOTES: 1. cr(x) CBM1 used known (x) 

2. a(x) CBM2 used = constant 

3. cf(x) C.F.E. is the conventional solution 
shown in Section III Alh 

4. a(x) C.Ave is the average at a node of 
the C.F.E solutions 



Error % is Exact-Approximate 

Exact 



x 100 
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Figure (12) shows three different approximations to the 
exact solution. The conjugate basis method with the exact 
stress' distribution gave the best solution at node 4 and 
had the least mean square error throughout. The next best 
approximation was the conjugate basis method where the 
exact stress distribution was not known. The conventional 
finite element solution is also plotted. 

3 . Linear Displacement Model and Constant Properties 

This is the same problem considered in the previous 
section only this time, k(x) = k - 1 . 

Equation (53) becomes 
ak L 

°N (x) = h dx *N (x) dx C56) 

The nodal stresses are 




The conjugate basis representation of the stress 

field is found using Equation (21) and is shown on Figure 

13. The exact solution in this case is a(x) = 2x. The 

conventional finite element solution is also shown. The 
ak 

factor - — 5 — has been disregarded. 
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a 

X 




Figure 13. o(x) vs x for k(x) = 1. 



4 . Uniform Strain 

Consider the same three element models with k(x) = 
600 to cancel the factors 1/40 and 1/15 resulting from the 
numerical calculations and U(x) = ^-g • This is a case 

where node 4 is displaced 0.3 units along the x axis. 



^ . V . ( x . ) = constant = .1 In this case the exact solution, 

dx * 

the conventional finite element solution, and the conjugate 
basis solution are exactly the same. This is shown in Fig- 
ure 14. This problem corresponds to a uniform bar under a 
concentrated end load P = 30. (EA) 
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Exact 

C.F.E. 

CBM 



30 



x 



a 



x 



Figure 14 • a(x) vs x for Uniform Strain. 

5 . Discontinuous Displacement 

In the previous examples, the stress was in fact 
continuous throughout the model. The CBM gave the best fit 
at the extremes and the least mean-square error over all. 
However, in problems where the stress is discontinuous due 
to a discontinuous displacement approximation, the CBM will 
not be the best fit. However, we consider examples with 
stress discontinuity to observe how the CBM tries to mini- 
mize the total error in a least squares sense in various 
models. The direct relation between discontinuous load and 
stress obtained here is true for the in-plane and bar prob- 
lems. However, other problems in mechanics may have other 
stress-load relations. For example, for bending structures 
(beams or plates), discontinuous stress arises from discon- 
tinuous moment, and not discontinuous force. 
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a. Displacement at Node 3 

For the same three-element bar given a displace- 
ment of 0.1 units at node 3, the displacement function is 
now 

U(x) = Jq 0 < x < 0 

fo 2 1 x < 3 (58) 

The strain for elements one and two is a con- 
stant 0.05 and for element three is zero. The exact stress 
distribution is given by cr(x) = k ( x ) where k(x) = 600. 

The conjugate nodal values are computed using 
Equation (23) and are shown below. 




600 

40 




(59) 



The conjugate stress field is computed using 
Equation (21) and is shown on Figure 15 along with the 
exact and conventional solutions. 
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Figure 15. a(x) vs x for Displacement at Node 3. 

The conjugate basis method spreads the error 
throughout the model to obtain the best least squares fit, 
thereby giving a stress field in the element where there is 
no stress. Consider the area under the stress curve as a 
measure of the average error. For the exact solution A = 

30 x 2 = 60 units; for the conjugate basis solution A = 

60.4 units. This gives a total percentage error of 0.66%. 
However, at node 3, the CBM value is only 53.3% of the maxi- 
mum stress at that point. Clearly the conventional finite 
element method gives a better approximation, 
b. Displacement at Node 2 

If the model had been displaced from node 2 
instead of node 3, the results would be similar in that the 
conventional constant strain element would be a better 
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approximation and the CBM spreads the error over the entire 
model. The results are shown below. 




Figure 16 . 



0 (x) vs x for Displacement at Node 2. 



Again comparing the areas under th 
curves A(exact) = 30; A(CBM) = 30.89. The err 
overall and 53.33% less than the maximum at no 

B. TWO-DIMENSIONAL IN-PLANE PROBLEMS 

The development of the computer program wh 
to solve the two-dimensional in-plane problems 
in Section IV. Appendix B gives a description 
ferent finite elements which were used to obta 
compare with the results obtained by the CBM. 
program itself is listed in Appendix C. 



e stress 
or is 2.96% 
de 2 . 
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The basic model used to verify the computer program was 
a 10 x 10 x 1 plane divided into either four or eight ele- 
ments as shown below. The units used for numerical calcula 
tions are inches and pounds. The figure also shows the 
global node numbers and the element numbering system. 




Model B 



Figure 17. Typical Two-dimensional Models. 

Each node has two degrees of freedom (DOF) , which are 
in-plane displacements u and v, and for a typical node n 
they are numbered as shown below. 



y 



DOF (2n) 



v = V, V 

n (x,y) 



n I DOF(2n-l)= u - u, N 

, n (x,y) 

i 

— i ■■ ■ 

n 
x 

Figure 18. Reference System. 

The following information is pertinent to all of the 
pr ob lems . 

E = 30 x 10^ lbs/in^ 
y = 0.3 

A 2 

All stresses plotted are G x (x,y) x 10 lbs/in 
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1 . 



Model A with Uniform Edge Displacement 




u(0,y) = 0. 

v ( 0 , y ) = 0. 

u ( 10 . , y ) = 1. 

v (10 . ,y) = 0 . 



Figure 19. a^CPSI x 10 ) for Problem III Bl. 



The plane is prevented from displacement in the v 
direction; therefore = 0, Equation (A-l) shows that the 
stress will be uniform across the model since the strain is 
uniform* This is the result that the CST solution gives. 
Because the CST solution is continuous for all elements, the 
conjugate basis solution cannot improve upon it. Both solu- 
tions in this case are exact for plane stress. 

This problem is analogous to the three element bar 
which was displaced at node 4. See Section III A4 . In the 
following figure, the stresses obtained by the conventional 
CST finite element solution are given in each element, while 
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the stresses obtained by the conjugate basis method (CBM) 
are the values given at the nodal points. 



2 . Model B with Uniform Edge Displacement 

This problem is analogous to both the problem of 
the bar pulled uniformly at its end (III A4) and the first 
Model A plane problem (III Bl) . It confirms that mesh A 
was sufficient for the problem and that the CST is a 
sophisticated enough element for this type of problem. The 
difference in boundary conditions between this problem and 
problem III Bl accounts for the difference in stress values 
In this case < 0, that is contraction in the v direction 
is allowed, thereby reducing 



300 



300 



300 
I 

I 




u ( 0 . , y ) = 0 

v ( 10 . , y ) = 1 



300 

I 



Figure 20. a (PSI x 10 ) for Problem III B2 
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3 . Model B with Uniform Edge Displacement, Modified 



Analytical values of 0^ were calculated for x 



xClOjy), x = x(5,y)and x = x(0,y)using the plane stress 
equations (A-l) . The values of strain were approximated by 
dividing the total contraction in the v direction for a 
given y by the original length 10. The strain ey was con- 
sidered constant along a constant y. The displacements 
were calculated in the CST routine. 



For x = 0 and x = 5 the CBM gives the better approxi- 



mation. At x = 10 the difference in percentage of error 
between the CST solution and the CBM is less than 0.67%. 

It is also interesting to note that the average of the CST 
values at nodes 2 or 8 is the same as the CBM value. This 
gives mathematical substance to the engineering practice 
of using the average value at a node. 



The values of a calculated were 

x 



a (0.) = 329 
x 



0 (5.) = 292 
x 



0 (10.) = 300 
x 



327 



309 . 



313 





311 



u(0 . , y ) 
v (0 . , y) 
u (10 . ,y) = 1 



327 



Figure 21. 
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4 . Model A with Linear Edge Displacement 



Here the right edge is given a linear displacement 

as shown below. Using the analytical procedure discussed 

4 

in Problem III B3, the values of a range from 329 x 10 

x 

along the edge (x,10) to zero along the edge (x,0). The 

CBM values are the best continuous approximation available 

for this model. A plot of 0 (0) vs y along (0,y) shows 

x 

how the CBM is better in approximating the maximum value. 



I 

I 

! u (0 . , y ) = 0 . 

' v (0 . , y) = 0 . 

I u (10 . , y ) = y/10 

I v (10 . , y) = 0 . 

I 

I 



24 

Figure 22. a (PSI x 10 4 ) for Problem III B4. 

; X 
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a Exact 
x 




Figure 23. O (0) vs y for Problem III B4 . 
x 

A comparison of the percentage of error at (x,y) = 

(0., 10.) that the CBM value is 7.4% too low while the CST 

average value is 28.5% low and the minimum CST value is 
50% low. 

5 . Model B with Linear Edge Displacement 

The' applied displacements in this problem make it 
similar to Problem III B4; however, the other boundary con- 
ditions in this problem allow for movement in the v direc- 
tion. In part the displacement solution of the CST program 
shows that at x=0., x = 5 . , and x = 10., < 0. Therefore, 

along the bottom edge where ex is very small, the effect of 

C becomes dominant and a negative a results, 
y x 
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From the results of the CST displacement solution, 
3 u 

the value of £ = t: — along the edge (x,0.) is 0.0026 and 

X O X 

3 v 

the value of £ = t: — = 0.016 along x=10. and 0.015 along 

y 3y 

(5,y) and (0,y). Using these values of £^ and the average 
of the e y' s Equation (A-l) for stress yields 

a (x,0) ' = - 73 x 10 3 PSI 

The CST solution is attempting to reflect this in 
elements 5 and 8. The CBM does show the compressive 
stresses and of course is continuous throughout the model. 
The average of the CST values at nodes where more than one 
element are incident is close to the CBM value at that 
node. For example, at node 5 the CST average is 150 and 
the CBM value is 148. This again confirms the existing 
practice of averaging. 



394 238 398 




u (o . , y) = 0 . 

u(10 . ,y) = y/10 



Figure 24. (PSI x 10 ) for Problem III B5 . 
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6 • Model B With Piecewise Linear Edge Displacement 



This problem has boundary conditions similar to 

problem III B 4; but, since it has two elements along 

x = x(10,y) and the applied displacement is only at node 1, 

the displacement profile of the edge is only piecewise 

linear with discontinuous slope. 

The expected stress profile is greatest at (x,10) 

where £ is maximum and compressive at (x,0) where £ =0 
x x 

but a compressive £^ exists. Likewise £^ is higher along 
(0,y) because £^ = 0 , than along (5,y). 

The CBM values conform with the expected distribu- 
tion and give a continuous stress field which is the best 
approximation in this case. 
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U (0 . , y) = 0 . 
v (0 . , y) = 0 . 
u(10.,0.) = 0. 
v (10 . , 0 . ) = 0 . 
u(l0. , 10 . ) =1 . 



Figure 25. 



O 



x 



(PSI x 10 



A 



for Problem III B6. 
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7 . Model A With Center Node Displacement 



For this problem the stress is discontinuous due to 
the displacement of node 5 in the interior of the model. 

The CST solution gives an exact representation of the stress 
field. Element 1 is in compression and element 3 is in 
tension . 

The CBM values are exact at the nodes. A stress 
plane constructed linearly thru the nodes would be in error 
at all ' locations other than the nodes. 

If a three dimensional plot is made with stress as the 
vertical axis, the volume under the CST plot would equal 
the volume under the CBM plot at 25 x O 

max 
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u (0 . , y) = 0 . 
v (0 . , y) = 0 . 
u(10. ,y) = 0 . 
v ( 10 . ,y) = 0 . 

u(5. ,5.) = 0.5 



Figure 26 . 



a (PSI x 10 ) for Problem III B7. 
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8 . Model A With Center Node Displacement, Modified 



The model is subject to a discontinuous stress field 
due to the internal displacement of the cneter node. If 
the model allowed more flexibility, one would expect the 
stress to be discontinuous in the region of node 5 and 
continuous along the edges. The CBM does not show the 
abrupt change at node 5 but does approximate the solution 
expected of a more flexible and hence more accurate model. 



329 -6 6 




u (0 . , y) 


o 

II 


v (0 . , y) 


o 

II 


u ( 5 . 5 . ) 


= 0.5 



Figure 27. a (PSI x 10 4 ) for Problem III B8 
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9 . Model B With Center Node Displacement 



This problem is similar to the previous problem on 
model A with the exception that the boundary conditions on 
this problem allow displacements in the v direction. The 
stress is discontinuous in the region around the center 
node, node 5. However it is unlikely that the stress is 
as severely discontinuous as the CST solution indicates in 
this region. Intuitively a has a maximum tensile value at 
(0,5) and a maximum compressive value at (10,5) and decreases 
in absolute value away from the line of action of the 
displacement . 

For this model, subject to a point discontinuity, 
the CBM gives a better overall representation of the stress 
field. 
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U (0 . ,y) = 0 . 

u(5 . ,5 .)= 1. 



Figure 28, 



a x (PSI x 10 ) for Problem III B9 
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10 . Model B With Uniform Centerline Displacement 



This is another problem with a discontinuous stress 
resulting from an internal displacement. Note that to 
minimize deflection in the v direction, the boundary condi- 
tion v(5,5) = 0 has been imposed. 

The expected value of a x (x) for x > 5 is zero. The 
CST solution is superior from two standpoints. It reflects 
the stress discontinuity very well arid it gives a better 
approximation to the stress field over the right half of 
the model. It is interesting to note that the average CST 
value at nodes 2 and 8 is the CBM value. 
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u (5 . ,y) = 1. 

u (0 . , y) = 0 . 
v(5. ,5.)= 0. 

-109 



321 154 -37 



321 154 




Figure 29. 



a (PSI x 10 4 ) 
x ' 



for Problem III BIO 
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11 • Model B Modified, With Uniform Centerline Displacement 



This mesh is a modification of model B and was 
designed to determine how rapidly the CBM stress would con- 
verge to zero. Each node along x=5 and x=10 was displaced 
an equal amount to force rigid body motion. 

The CST model responds to the discontinuity quickly 
as expected. The CBM does try to converge and is a better 
approximation over this part of the model than the approxi- 
mation with the previous mesh (III. B. 10). The conclusion 
is that if a discontinuity occurs across an entire model, 
and there is sufficient room for the continuous function to 
correct itself, and the area of interest is not at the 
discontinuity, the CBM will be the best approximation for 
a continuous field. 




U (10 . , y) = 1 . 
u (5 . , y ) = 1 . 
u (0 . ,y) = 0 . 



Figure 30. a x (PSI x 10 4 ) for Problem III Bll 
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12 . Model B Tapered With Uniform Edge Displacement 



This model is another modification to model B. The 
length has been doubled and a taper introduced by making the 
right edge one half the height of the left edge. Symmetry 
about the horizontal centerline has been maintained. The 
model simulates a tapered plane and is the two dimensional 
analog of the one dimensional tapered bar. Both the CST 
solution and the CBM give the expected stress distribution. 
However the CBM is continuous over the entire model and 
therefore is a better approximation. 




u ( 2 0 . , y ) = 1. 
U (0 . ,y) = 0 . 
v (0 . ,y) = 0 . 



Figure 31. 



O (PSI x 10 4 ) for Problem III B12 
x 
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13 . Composite Material Models With Uniform 

Edge Displacement 

This problem consists of using the CBM with a CST 
input on the composite material models used by Lin, Salinas, 
and Ito [Ref. 6] and comparing the results for stress 
convergence and computer time usage. The models used in 
Ref. 6 are shown in figures (32) and (33). The PLISOP model 
shown in figure 34 has been used by Salinas subsequent to 
publishing Ref. 6. The stresses resulting from using the 
CST and LST routines on figures (32) and (33) and the PLISOP 
routine in figure (34) (see Appendix B) are shown on figures 
(35, (36), and (37) along with the CBM values. 

In figure (37), along the edge (5.2,y), where there 
is only one material present, the CBM yields a smooth curve. 
However, where two materials are present and there is a 
sharp change in stress due to the different Young’s Modulus 
of the two materials, the CBM tries to smooth out the change 
in a manner similar to the discontinuous problems shown 
earlier. As a result several CBM model values in the region 
of the change in material are not in agreement with values 
obtained from" the other solutions. The results are in good 
agreement away from this area of change. 

An analysis of the amount of computer time each 
program used is shown in Table II. 
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MODEL 3 MODEL 4 LST PLISOP 

(124 DOF) (242 DOF) MODEL 3 16 ELEMENTS 

(430 DOF) (204 DOF) 
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Figure 32.' Finite Element Grid System, 124 DOF. 
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Figure 34. Finite Element Grid System for PLISOP, 204 DOF. 
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a (x 10 5 ) 

X 




Figure 36. a (x=y) (PSI x 10 5 ) for Problem III B13. 

X 
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a (x) (x 10^) 
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IV. THE COMPUTER PROGRAM DEVELOPMENT 



In developing a computer program to determine stresses 
utilizing the conjugate basis method, four major requirements 
had to be met: 

i. The inverse of the global fundamental matrix is needed. 
This requires forming the element fundamental matrices 
and merging them into the global fundamental matrix 
itself . 

ii. An' approximation to the stress field in the element is 
required . 

iii. The conjugate nodal stress values for each element are 
formed . 

iv. The conjugate basis approximation to the stress field 
must be assembled. 

There are many existing programs which satisfy the 
second requirement. In this thesis a program utilizing the 
CST written for use in Ref. 6 has been used. It was chosen 
because it was available; it uses a simple two-dimensional 
element — the three noded triangle; and problem solutions 
using this element were available for comparison as to the 
computer time used and the accuracy of results. It was not 
an objective of the thesis to develop an optimum conventional 
finite element solution. Any other CST program or a more 
sophisticated program could have been used. The advantage 
of the CBM is that it makes the stress at a node a function 
of the entire model and not just the element. 

Therefore the development of the CST program itself will 
not be discussed with the exception of a brief comment on 
the constitutive relation matrix, C. For plane stress 



70 



problems where 0 = T = T = 0 

z xz yz 




1 u 0 

y 1 0 

0 0 . 5 (1-y) 



where {a} = [C]{e} 



For the CST, each displacement field is approximated by 
a linear field over each element, hence the local basis 
functions are area coordinates. The local fundamental 
matrix, which is the inner product of each of the basis 
functions, is therefore dependent on the total area of each 
element. The integration can be found in Felippa [Ref. 7] 
and the resulting matrix is 

2 11 
12 1 
112 

After the area of each element is calculated in the 
normal CST routine, the element fundamental matrix — EFM — 
is formed. It is immediately merged into the global funda- 
mental matrix — GFM. The merging routine is described in 
section III.A.l.d. 

After the entire GFM is formed, it is inverted using 
subroutine SYINV. This is the most time consuming addition 
to the CST program. See Table II. Notice that for the 92 
element model SYINV used 20 seconds but for the 200 element 
model the time was 133 seconds. 



A 

12 



71 



The conjugate basis nodal stresses are calculated in 
two steps at the completion of the standard CST program. 

The assumption is made that stresses will only be calculate 
at nodes. Since <J>^ « 1 at a node (see section II. C. and 
section III.A.l.e.) the values of the local conjugate basis 
function at a node are just elements of the GFM inverse, 
in accordance with equation (44). This assumption is not a 
serious limitation because the finite element mesh can be 
arranged to have a node at any point of interest or the 
stress can be calculated by interpolating between neighboring 
nodes. However it greatly simplifies the computer 
calculations . 

Recall that the conjugate nodal values of stress for 
each element are given by Eq. (23) as 

L 



The CST value of a(x) is constant for each 

taken from the CST solution. The solution 

A ( e ^ 

is found in Ref. 7 and is — - — for each N. 
conjugate value for each node of element e 



element and is 
L 

of / 



^ N (x)dx 



The ref ore 
A ( e ) 



xs 



the 

0(e) . 




A (e) 
3 



0(e) 



(60) 



The conjugate function representation of stress is given 
by Eq. (21) as 

- l oj 

e = l 
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In the computer program a(x) is calculated one node at 
a time. The value of 0^ is multiplied by the element of the 
GFM inverse in the row corresponding to the global node that 
the local node N is incident to and the column corresponding 
to the global node being calculated. This multiplication 
is carried out for each element and the resulting sum is 
printed out as the nodal stress. 

The modification of an existing conventional finite 
$ 

element program to perform conjugate basis calculation of 
stress requires the. addition of three simple algorithms, 

i. An algorithm which forms the global fundamental 

matrix (GFM). This algorithm may be located in the 
element stiffness DO loop, where the element funda- 
mental matrices (EFM) are individually formed and 
whose elements are inserted in the appropriate place 
in the global fundamental matrix, C^p. 

ii. An algorithm to invert the global fundamental matrix 

AT 

C^p and obtain C (see equation (12). C^p is a banded 
symmetric matrix and inversion algorithms are available. 
(In this work E. Wilson's SYINV subroutine was used.) 
This algorithm is located after C^p is formed, and is 
the most time consuming aspect of the conjugate basis 
method. 

iii. An algorithm forthe calculation of nodal stresses 
associated with equations of the conjugate basis 
method. This algorithm uses the conventional stresses 
of the finite element method and must therefore follow 



the conventional stress calculation routine. 



A listing of a CST finite element routine is given in 
Appendix C with the CBM modifications shown enclosed in 
b oxes . 
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V. CONCLUSIONS AND RECOMMENDATIONS 



As a result of this study several conclusions can be 
made concerning the use of conjugate basis functions with 
in- plane problems . 

Where the model is not subject to interval discontinu- 
ities, the CBM gives the most accurate approximation to the 
stress over the entire model for a given displacement. 

If there is a point or localized discontinuity, the CBM 
will be a poor representation of the stress near the dis- 
turbance but will be a good approximation away from the 
disturbance. Mesh refinement can localize the disturbed 
area. It is understood that points of stress discontinuity 
are usually points of interest and that mesh refinements 
generally increase computer time, therefore the use of the 
CBM is limited. 

If there is a stress discontinuity across the entire 
model, the CBM is a poor choice if the model is too small 
or too coarse to allow the CBM to settle out. Refining the 
mesh may be worthwhile if a continuous approximation to the 
stress is wanted in regions away from the discontinuity. 
However the additional computer time needed places a 
premium on this marginal gain. 

If a large number of problems are run on the same model 
and the global fundamental matrix can be inverted once and 
stored, the additional computer time per problem for the 
CBM solution is minimal. The advantage of a continuous 
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stress solution can be had without using the more sophisti- 
cated elements with nodal strain degrees of freedom. 

Finally, at nodes where the stress should be continuous 
and several elements were incident at that node, the average 
of the CST stress values at that node agreed very closely 
with the CBM stress value. This gives mathematical sub- 
stantiation to an engineering practice. 

It is recommended that 

i 

1. The CBM be used for in plane problems without interval 
discontinuities because the resulting stress field will 
be continuous, which the actual stress field is, and the 
stresses will be consistent with the displacements. 

2. Further studies with the CBM should be carried out on 
other classes of problems. For example, plate problems 
or beam problems might be better approximated with CBM 
solutions . 

3. Additional studies be carried out to determine the cost 
effectiveness of incorporating and using the conjugate 
basis method in existing conventional finite element 
programs . 

4. An investigation be undertaken to determine if the CBM 
can be modified to yield stress discontinuity when 
appropriate . 

5. A study be made to establish the relation of domain 
modeling to CBM analysis. 
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APPENDIX A 



THE IN-PLANE PROBLEM 



In this discussion a plane is defined as a body bounded 
by two surfaces of zero curvature whose thickness, h, is 
mush less than its surface dimensions. For in-plane 
problems only axial loads N(x,y) are considered. There are 
no lateral loads or moments. Displacement boundary condi- 

i 

tions can be of the type au 4- bv = c where a, b, and c 
are specified coefficients; and u and v are displacements. 
In plane forces may also be specified on portions of the 
boundary where displacement is not specified. 




Since k is small in comparison with the plane dimensions, 

O << O >0 and x , x < < t . This suggests the use of 

z x* y zx zy xy 00 

plane stress conditions, that is, a , T , T are considered 

z z x z y 

to be zero. 

The basic stress strain relations become the plane stress 
equations . 
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o = r (e + y e ) 

x , 2 x y 

i - y 



a = -z (ye + e ) 

y i-p 2 x y 



T 

xy 



E 

1 + y 



G 

xy 



The strain displacement relations are 



G 

x 



9 u 
9 x 



G 

y 



9 v 

9y 



G 

z 



, 9 u 

( 37 + 



9 v 

9y 



) 



1 e 3u 9v N 

G ** 7T ( ~ — + T — ) 
xy 2 9y 9x 



Substituting Eqns • (A-2) in Eqns. (A-l) , the partial 



differential equations for the-in plane problem become 



a 

x 





u 



9 v 

9 y 



] 



a 

y 




9 u , 9 v -i 

u 37 + 37 ] 



T 



xy 



E r 9u 9jv -I 

2(1 + y) L 9y 9x J 



(A-l) 



(A-2) 



(A- 3) 
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Stress resultants N , N , N having the units of force 

x y xy ° 

2 

per unit length (// in/in ) are defined as 



h/2 

N = f adz 
x J x 

-h/2 



N 



xy 



h/2 

/ 

-h/2 



T d z 
xy 



N 

y 



h/2 

/ a y dz 

-h/2 



(A-4) 



These last expressions imply internal discontinuity in 
stress is associated with a discontinuity in applied internal 
in plane-force. 

For plane stress the equilibrium equations are 



I F 



x 



0 




3t 



xy 

3y 







3N 



xy 

3y 



0 



(A-5) 



1 F 



y 



0 : 




(A-6) 



The finite element method yields a system of algebraic 
equations which are obtained from the above partial 
differential equations by the usual method as given by 
Zienkiewicz [Ref, 8], 
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APPENDIX B 



THE FINITE ELEMENT METHOD 

An extensive description of the finite element method 
can be found in Zienkiewicz [Ref, 8]. However to facilitate 
understanding of this project, certain features of the 
elements will be presented. 

In all finite element programs certain rules must be 
followed in selecting the assumed displacement field. 

1. The finite element displacement field must accommo- 
date rigid body displacements. 

2. At common boundaries two elements remain in 
continuous contact with each other. Neither gaps nor over- 
lappings are permitted. 

3. Constant straining patterns must be included in 
higher order descriptions of displacement. 

CONSTANT STRAIN TRIANGLE (CST) 

For the CST the u and v in-plane displacements are 
approximated by a linear field over each element. Since 
strain is the first derivative of displacement, this gives 
a constant strain element*. 
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The area coordinates used in the program are defined 



as 




With these coordinates the linear displacement field 
over the element is simply 

u(x,y) = L x + L 2 u 2 + L 3 u 3 

v(x,y) = L x v x + L 2 v 2 + L 3 v 3 (B-l) 
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r 





LINEAR STRAIN TRIANGLE (LST) 



3 




Figure B- 2 . Linear Strain Triangle 



The LST. has three nodes on each side and each displace- 
ment field is taken as a 2nd order polynomial over the 
element domain. The displacements u and v for each element 
are as follows 

u(x,y) = N 1 u 1 + N 2 u 2 + N3II3 + N^u^ + N 5 u 5 + N 6 U 6 

(B- 2 ) 

v(x,y) = U 1 v 1 + N 2 v 2 + N 3 v 3 + + N 5 v 5 + NgVg 

where = 2L^ - 1 for a typical corner node 

= 4 L^L 2 for a typical mid side node 

The are called interpolation functions for the 
assumed displacement field. Differentiating the displacement 
fields gives linear strain over each element, hence the name 
linear strain triangle (LST). 
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PLISOP 



The PLISOP program was coded by Professor Gilles Cantin 
and his students at the Naval Postgraduate School, Monterey. 
It performs a plane stress or a plane strain analysis by 
the finite element method using numerically integrated 
isoparametric elements. These elements have four sides, 
each of which may be either curved or straight. The same 
shape functions which define the geometry also are used to 
describe the displacement. The strains will be a function 
of the number of nodes per side of the element. The 
comparisons made in this thesis were to a cubic element. 

This isoparametric element uses a cubic displacement field 
and yields parabolic strain (stress) fields over each 
element . 



All of these elements, the CST, LST and PLISOP, with 
displacements being the only degrees of f reedom, yield 
discontinuous stresses along element interfaces. The CBM 
gives continuous stress along element boundaries. 




Figure B-3. 12 Node PLISOP Element 
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For each of the conventional finite elements, the order 
of the strain approximation is one less than the order of 
the displacement approximation. However for the CBM the 
displacement approximation and the strain approximation are 
of the same order. 
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APPENDIX C 
COMPUTER PROGRAM 
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READ (5,30) ( ( I T E M ( I » J ) , J = I , 3 ) ,T(I),I = I,NUMEL) 
WRITE* 6, 35) ( I , T ( I ) ,(ITEM(I,J),J = 1,3) ,1=1, NUMEL) 
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